function hydro_run()

  jobname = 'turb1-devel';
  t_end = 10;  % ref 1
  Nx = 32;
  L = 1500; % ref 1
  Temp = 1; % ref 100
  hydroOrder = 2;


  t_steps = [0:1:10];
  params = hydro_params(Temp, Nx, L, t_steps, hydroOrder);

  params.plotHydro = true;

  params.RelTol = 1e-8;
  params.AbsTol = 1e-8;


  % output filename with current timestamp
  datestring = datestr(clock, 'yy-mm-dd-HHMM');
  fname = sprintf('%s-%s-T%d-L%d-Nx%d-tbeg%d-tend%d', datestring, jobname, round(Temp), round(L), Nx, 0, round(t_end));

  % initializes derivative kernels, plotting 
  hydro_init(params);

  % initial conditions, or replace with your own
  uIn = hydro_ic(params);

  % run the ode solver
  [tOut, uOut] = hydro_solve(params, uIn);

  % save results
  save(fullfile(pwd, fname), 'tOut', 'uOut', 'params', '-v7.3');
end
